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Abstract. Any eigenfunction of the laplacian on the sphere is given in terms of a 
unique set of directions: these are Maxwell's multipoles, their existence and uniqueness 
being known as Sylvester's theorem. Here, the theorem is proved by realising the 
multipoles are pairs of opposite vectors in Majorana's sphere representation of quantum 
spins. The proof involves the physicist's standard tools of quantum angular momentum 
algebra, integral kernels, and gaussian integration. Various other proofs are compared, 
including an alternative using the calculus of spacetime spinors. 



1. Introduction 

When manipulating functions on the sphere, it is convenient to have a functional 
representation whose behaviour under spatial rotations is transparent. This is not 
provided by the spherical harmonics Yp{6, (p), which are required to transform as basis 
functions for the appropriate matrix representation of the rotation group SO (3). In his 
Treatise on Electricity and Magnetism |Tj, chapter IX, Maxwell found an approach to 
real sums of spherical harmonics, representing such functions as directional derivatives 
("multipoles") of 1/r; the resulting homogeneous polynomial in x,y,z is the function 
represented by spherical harmonics on replacing the cartesian coordinates with the 
appropriate polar expression. The directions of these derivatives rotate directly as 
vector-like objects. Sylvester |2] proved that Maxwell's multipole representation is 
unique, this result (amongst others) being known as Sylvester's theorem. The theorem 
has been discussed and proved by several authors: Courant and Hilbert [Sj using 
Bezout's theorem from algebraic geometry; Backus ^ using the first isomorphism 
theorem for rings; and Zheng and Zou [5] using traceless symmetric tensors. 

In this paper, I aim to provide an alternative, concrete proof of this result, using 
mathematical tools in every physicist's kit |: quantum angular momentum algebra, 
integral kernels, and gaussian integration. In section [3 comparison is made between 
these various proofs. The main object utilised is Majorana's sphere representation for 

I It is unlikely that Sylvester would have approved, given his disparaging comments about Maxwell's 
method at the end of his note. 
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quantum spins [Uj. In this representation, a spin state (of definite spin j) is represented 
by 2j points (unit vectors) on tfie Riemann sphere, generahsing the famihar Bloch sphere 
representation for spin 1/2. Majorana himself introduced his sphere representation to aid 
calculation of transition probabilities in quantum mechanics jH] • It has been popularised 
more recently by Penrose [Zj, [Hj (particularly appendix C), and has been used to describe 
some Bell-type inequahties to quantify chaos and calculate geometric phases for high 
quantum spins piU ITT] , and characterise polarization in 3-dimensional electromagnetic 
and other fields ^21 The main objective of the present work is to demonstrate 
that in the special case of a real state of fixed spin, the Majorana vectors are the 
Maxwell multipole directions. Majorana's representation is more general, applying to 
both integer and noninteger spins and complex functions, but does not, in general, have 
the concrete multipole interpretation described below. 

Maxwell's multipole representation represents a real function f{6,(j)), with definite 
angular momentum + - that is, it is eigenfunction of the laplacian on the unit sphere 
with eigenvalue — j(j + 1) - as the direction-dependent part of a multipole derivative; 
i.e. there are j unit vectors Ui, . . . ,Uj such that 

f = D^,...D^^^, (1.1) 



where D^. represents the directional derivative operator Ui ■ V, and r = ^y x"^ + y'^ + z"^ . 
It is, of course, straightforward to show that a function of this form is an eigenfunction 
of the laplacian (see, for instance, |14J, section 28); however, it requires more work to 
show that any eigenfunction can be put into this form - this is Sylvester's theorem. 

A more useful form of (jl.ll) is obtained by finding the derivatives explicitly, then 
expressing geometrically in terms of vectors (e.g. via a homogeneous polynomial 
representation in cartesian coordinates x,y,z). In effect, the directional derivatives are 
replaced by scalar products with the corresponding vector, and there is an additional 
term consisting of times a complicated combination F of the components of the Ui, 

f{9, 4>) = C{u{9, 0) • ni) . . . {u{9, 0) ■ uj) + r^F. (1.2) 

Here, the unit vector u{6, 0) is the point on the sphere given by the spherical coordinates 
6, 0, and here and hereafter, C is an unspecified constant factor. The equivalence 
between the representations (jl.ll) and p.2|) is demonstrated in Appendix A| using Fourier 
integration methods. The latter form of Sylvester's theorem will be proved here. 

The Maxwell multipoles have simple forms for the real parts of the spherical 
harmonics YJ^{6,(j)): for zonal harmonics (m = 0), the j multipoles are aligned in 
the z-direction, whereas the tesseral harmonics (m > 0) have m multipoles in the xy- 
plane arranged in a regular 2m-gon, and the remaining j — m aligned in the z-direction. 
A further simple example is the case when j = 2, where the function is equivalent to 
a traceless symmetric tensor. On adding a sufficiently large, arbitrary constant times 
the identity matrix, / can be represented geometrically by an ellipsoid: the two planes 
in which the projections of the ellipsoid are circular are perpendicular to the Maxwell 
multipoles in this case. 
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2. Quantum notation, spherical harmonics and spin coherent states 

Notation from the quantum mechanics of spin will be used, although it applies for 
general functions on the sphere. This provides a convenient representation of the linear 
algebra of the underlying mathematical structures, i.e. the irreducible representations of 
the rotation group S0(3) and the special linear group SL(2,C). {ip) denotes an abstract 
ket state with definite integer or half-integer spin j, fixed once and for all (i.e. it is 
an eigenket of total angular momentum J^). j will be assumed to be an integer unless 
stated otherwise. 

With respect to a certain coordinate frame, may be written 
j 

l^)= Yl (2.1) 
m=-j 

where \m,j) is an eigenket of the angular momentum component operator with 
eigenvalue m, and where no confusion will ensue, \m,j) is written |m), and the same 
applies for the corresponding bras. The phase convention used is the standard Condon- 
Shortley convention given by applying successive raising operators to | — j). These 
states will often be rotated by 3-dimensional rotation operators; the rotation operator 
taking the unit vector in the +2;-direction, u{0), to the unit vector u{6,(j)), is denoted 
Rz{4>)Ry{6). In terms of Euler angles, the vector is first rotated through 6 in the y- 
direction, then through in the z-direction. The notation will be abused, and rotation 
of the corresponding bra-ket states will be represented by the same operators: 

R,{(l))Ryie)\m) = \m;e,(l)) (2.2) 

where the ket on the right hand side is an eigenket of the angular momentum component 
operator u{6, 0) ■ J with eigenvalue m. 

It may be shown, either from group theory ^Hj, or other methods fHlE!, that for 
j integer or half integer, the rotation operator matrix elements (Wigner D-functions) 
are given by 

= exp(— im^)^/ (j + m')\{j — m')!(j + m)!(j — m)\ 

/ 1 \k—ra'+m„2j+m'—m—2k„m—m'+2k 

X ^ ' - f (2.3) 

^-^ [m — m' + k)\(j + m' — k)\(j — m — k)\k\ 

where in the sum, c = cos ^/2, s = sin ^^/2 and the index k runs through all integers with 
nonnegative arguments in the denominator factorials. Functions of spatial coordinates 
(i.e. in 'coordinate representation') depend on these matrix elements, as we shall see. 
Here and in the following, the active rotation convention of ^3 will be used, rather 
than the passive convention of [13 CHI UH] • 

A standard result (see, for example, [IH] equation 4.17a, |T2j equation (3.6.52)), on 
putting m' = in ()2.3|) and conjugating, gives the 'fully folded' representation of the 
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ket |m) as a spherical harmonic, 



47r 



(0; 9, 0|m) = V%{<peor = ^ ^f^Ypi^^, <P)- (2-4) 

The term 'fully folded' refers to the fact that the m' = state has been used in this 
representation; it is equivalent to mapping the spin states into basis functions for the 
rotation group SO (3). This result gives an equivalence between general spin j states 
and functions decomposed in terms of spherical harmonics: 

47r 

^— Y (2-5) 

m=~j 

Directions on the sphere, instead of being represented by the spherical coordinates 
9,(j), can be represented by the stereographic coordinate in the complex plane 

C(^, 0) = exp(i0) tan^/2, (2.6) 

which is a complex number (possibly oo) obtained by stereographic projection from the 
south pole. When the context is clear, this coordinate will be denoted simply by C- It 
will often be useful to use this stereographic representation; when appropriate, |m; 9, (p) 
shall be written |m; (), etc. It is easy to find ()2.3|1 when m' = — j (and j is integer or 
half-integer), giving the spin coherent state |2()j : 

C _ I y+m^j-m^j+m 

(-j;C|m) = exp(imargC)\/(2j)!(j + m)\{j - m)\—— -— — 

[J + m)l[j — my. 

c^' exp(-ij arg C) (-1)^'+'" [exp(i arg C)C]''+" 



(j + m)\{j — m)\ 
exp(-ijargC) 



:i + ici 



2V 




Thus, in the 'fully extended representation' (choosing m' = —j), \m) is represented as 
an m-dependent monomial of ( with positive power, times an m-independent function 
of ( and a numerical factor. For brevity, this numerical factor is written 

1/2 

^ (-1)^+™ I 1 . (2.8) 




In the fully extended representation, states are written as coherent basis functions for 
the linear group SL(2,C). A general ket is represented by its Majorana function [HllHl 

p40 = {-r,m 

^ ISI / m=—j 

The function p^{C) is a polynomial in (, times a prefactor which depends nonanalytically 
on (, and all |?/')-dependence is in the polynomial part. The polynomial will be called 
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the Majorana polynomial; the zeros of p^(C) are the zeros of the polynomial. If m' were 
taken to be +j rather than —j in ()2.7|) . the polynomial in ()2.9p would depend on (*, not 
(; this representation was used by Penrose with north-pole stereographic projection. 
Unlike the fully folded representation ()2.5p . the Majorana representation can be used 
when j is an integer or a half integer. 

The abstract ket {ip) can therefore be represented by the fully folded state ()2.5|) . 
given in terms of spherical harmonics, or the fully extended state, whose details depend 
on the Majorana polynomial in the stereographic coordinate (. It is useful to transform 
between these two representations, and to do this, we introduce the extending kernel 
(where ( is assumed independent of 9,(j)), 

K{c\e,<p)^'^^{-r,C\0;e,<p) 



2j + l 
Air 



J2 {-j;C\m){m\0;e,(P) 



^^ITTW^ ^ vr(M).„C-. (2.10) 

V IS I / rn=—j 

In the second line, a resolution of the identity has been used, and in the third, ()2.5p 
and ()2.9|) . The integral kernel can therefore be written as a Majorana function, 
whose coefficients are the conjugate spherical harmonics, depending on 6 and 0. It is 
straightforward to verify that K{(\6, 0) indeed transforms between the representations: 

rd<p r de smei;{e,<j))K{c\e,<f)) = p40. (2.11) 

JO Jo 

Its inverse, the folding kernel K'^{6,(j)\() is defined 

K\e,<j>\o = ^^{0;e,<j>\-r,O- (2.12) 

This kernel provides the inverse of ()2.11|) : 

,,^^^^p4C)K\e, 010 = ^{e, <p). (2.13) 



3. The Majorana sphere 

The equation for the Majorana polynomial ()2.9p applies for any integer or half-integer 
j. In particular, if j = 1/2, the Majorana polynomial is linear in (. Its single root is 
?/'_i/2/V'+i/27 which may be inverse stereographically projected onto the Riemann sphere. 
In quantum mechanics, this representation of a spin 1/2 spinor by a point on the unit 
sphere is known as the Bloch sphere representation; it is useful because the point on the 
Bloch sphere is independent of the overall normalisation and phase of the state. 

The Majorana sphere is the generalisation of the Bloch sphere construction to 
polynomials of arbitrary order 2j. By the fundamental theorem of algebra, the Majorana 
polynomial has 2j complex roots, which are inverse stereographically projected onto the 
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Riemann sphere, giving 2j unit vectors, which are points on the Majorana sphere. These 
points rotate rigidly when the state is rotated, and such rotations correspond to unitary 
Mobius transformations on ( in the Majorana function. These unit vectors will be called 
Majorana vectors. 
Therefore, 

MO - t ■M*'" 

^ ' ' ' m=—j 

exp(-ijargC) ^ .2j , ( ^ /- \ a ^\ 

= (l + |^|2), (-1) 11(C - Cn), (3.1) 

where the 2j roots are labelled Cn, with spherical coordinates 6'„, 0„. The terms in this 
product are all linear in (, and may therefore be viewed as spin 1/2 matrix elements. 
Including the appropriate normalisation, and writing in the angular momentum j, this 
means that 

^ ^ :T7in(-2'2'^ "2' 2'^")- ^^-^^ 

n=l ^ ' 

This may therefore be called Majorana's theorem. It states that the coherent state 
wavefunction of a given state of spin j (integer or half integer) is proportional to the 
product of 2j spin 1/2 coherent states. The proportionality factor depends on ipj (whose 
modulus may be expressed in terms of the symmetric polynomials of the Majorana 
polynomial, and their conjugates), and the roots C„. It is invariant with respect to 
unitary Mobius transformations of the roots, thus the right hand side of ()3.2|) rotates 
rigidly, according to the rotation of spin 1/2 states only. This result does not appear to 
be particularly well-known in the literature of quantum angular momentum. 

An important observation, made by |3J |3] (for polynomials analogous to the 
Majorana polynomial), is that the Majorana vectors form antipodal pairs for a state 
whose fully folded representation is real, as will now be shown. Let f{6, (p) be a real 
function on the sphere, with definite /, given in terms of spherical harmonics 



m=-j 

Spherical harmonics are related to their conjugates by the following standard expression: 



Yp{e,ct>) = {~irY-^{e,cl^). (3.4) 

Reality of / therefore requires that the complex coefficients satisfy an equivalent 
condition: 

al, = {-lTa_^. (3.5) 

The Majorana function p/(C) corresponding to / is given by ()2.9j) . with ipm replaced 
by am. The stereographic coordinate representing the antipodal point to C is — 1/C*, as 
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( , /^*\ exp(-ijargC)|Cp-^' . ,v+m^*-j 



m=-j 
j 



may be directly verified from (j2.6|) . Putting tliis into tlie Majorana function, 

= '""^Vl+lcpi^^^' ^ (-l)"^ar./i™|CP-'"exp(i(j +m)argC) 

V ISI j m=—j 

_ exp(ijargC)(-lV >^ ^ 

V ' ' ' rn=—j 

= i-iypficr, (3.6) 

where in the third hne, the identity ()3.5|) and the fact that = fim have been used, 
and in the fourth, — m has replaced m in the summing index. This result shows that ( is 
a root of pf if and only if — 1/C* is; the Majorana vectors of real functions are antipodal. 
Since the extending kernel K{(\6,(j)) can be written in Majorana form ()2.10p . with 
coefficients Yp satisfying (j3.5p . this shows that the Majorana vectors of the kernel are 
arranged in antipodal pairs. 

In quantum mechanics, such 'real' states occur in the presence of time reversal 
symmetry, since it can be shown that the time reversal operator antipodises the 
Majorana vectors of the state |13j. Thus, in situations such as Kramers' degeneracy 
|19j . when j is an integer, the states are real and there is no degeneracy, whereas if j 
is a half integer, the states are twofold degenerate with mutually antipodal Majorana 
vectors. 



4. The Majorana polynomial for spin 1 and 3-dimensional vectors 



The simplest Majorana function which has a folded representation occurs in the case 
j = 1, when the polynomial is quadratic. In this case, 

exp(— i argC) 



Pi'iC) 



(1 + ICPV 

exp(-iargC)V'+i 



(C-C-)(C-C+) 



(4.1) 



(1 + ICP)^- 

where the roots (± are 

> ^ V^o ± V^o - 2V'+iV'-i 

This can be seen more clearly if the state represented by a spherical basis \m) is 
replaced by a (generally complex) vector v in a cartesian basis, according to the standard 
transformation 



(4.2) 





( 

1 


Vy 




\v. J 



-i 







^+1 


\ 






^0 




/ 


v 


^-1 


/ 



(4.3) 
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From here, it is clear that v is real if and only if = and V'o is real. Since the 

matrix is unitary, v is normalised: ■ t> = 1. In a cartesian basis, the pair of roots ()4.2p 
becomes 



a = (4.4) 

-Vx + Wy 

If the state is real, with spherical coefficients a+i,ao,a_i satisfying ()3.5|) (so v is 
also real), the expression for pair of roots simplifies: 

^ V2a-i ^ + iVy 

This is the familiar stereographic representation for the pair of antipodal unit vectors 
±.v. This imphes that any real state of spin 1 corresponds to a ket |0, 1; 6, 0), where 6, (j) 
are the spherical coordinates of one of the pair of antipodal Majorana vectors; this is a 
real, 3-dimensional vector up to a sign. 

It is therefore natural to ask what the ket | + 1, 1; corresponds to in 3- 

dimensional geometry. Such a state is a rotation of the basis state | + 1, 1), which, by 
fl2.7|l . has Majorana function 

(-M;CI + i.i)^^^5^(ziHiOc=. (4.6, 

This has a repeated root at zero, corresponding to a pair of Majorana vectors in the 
+2;-direction (thus these spin representations can be used to define spin 1/2 spinors 
|21j). Therefore, the two Majorana vectors for | + 1,1;C(^,0)) are the same, and are 
in the direction u{9, 0) since the Majorana sphere rotates rigidly. The case of repeated 
roots corresponds to the discriminant of the Majorana polynomial being zero; by ()4.2p 
and ()4.4j) . this implies that, in a cartesian frame, v{6,(j)) ■ v{6,(f)) = (note there is 
no conjugation). Such vectors are called nilpotent (or isotropic) fn\ I22j: the nilpotent 
cartesian vector corresponding to | + 1, 1; () will be denoted v{C)- K can also be shown 
1221131 that the real vector x is parallel to the ^-direction; such vectors in 

the electromagnetic field E + iB have relativistically distinguished status (e.g. |23l EH] 
and 123 p258). 

Any inner product involving spin 1 states can therefore be written in terms of inner 
products of cartesian vectors (which may be real, nilpotent or general complex). This 
implies that the vector corresponding to the ket | — 1, 1; () is ^'(C)*; giving the required 
identity (-1, 1; C| + 1, 1; C) = ^^(C) " ^^(C) = 0. 

For any real state, corresponding to \0,l;6,(f)), the normalisation + |aoP + 

|a_ip = 1 implies that, up to a phase factor 

a^, = ^Y,'*{e,<py, (4.7) 

this constant appears as a multiplying factor in the Majorana function ()4.H) . and will 
be used in the next section. 
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5. Analogue of Majorana's theorem for real states 

The form of the Majorana polynomial for general basis states is given by ()2.7j) : the 
polynomial for \m) has j + m vectors in the +z direction n(0), and j — m in the —z- 
direction ^(Tr). As described in the previous section, since the Majorana sphere rotates 
rigidly, the Majorana function |m; 6, 0) has j + m Majorana vectors in the direction 
u{6, 0), and j — m vectors in the opposite direction —u{6, 0) = u{7i — 6,(f) + -k). When 
j is an integer, the extending kernel K{(\6,(j)) of ()2.1()|1 exists, and by its definition as 
a Majorana function, it has j repeated Majorana vectors in u{6,(j)) and j in —u{6,(j)). 
Therefore 

4^ r^fnp ^\ - exp(-ijargC) v"^*(ft 
exp(— ij arg^)Y'^*(6' '^'^ 



;i + ici 



' ■(-1,1;C|0,1;^, 



{V^Yi*{e,(i,)/2^y 

= C(n(^,0)-KC)y. (5.1) 
As before, C denotes an unspecified numerical constant, and inner products have been 
rewritten in a cartesian basis. In the penultimate line, the appropriate normalisation 
factor for the j = 1 polynomial from ()4.7|) was used, and the 6', dependence in the 
prefactor here cancels. A consequence of the factorisation in the second line implies 
that the conjugated spherical harmonics Y^* may be constructed as Y-* times the 
appropriate symmetric polynomials in the roots C,{9, 0) and — 0)* and a numerical 
factor. The final line of shows that any extending kernel for integer j is proportional 
to the jth power of the dot product of a real vector and a nilpotent one. 

The Majorana polynomial p/(C) for the real function / fl3.6|l . may be written in 
a similar way to (jS.lj) . since its Majorana vectors form antipodal pairs. Choosing and 
labelling a vector Un = u{6n, 4>n) from each pair, with n running from 1 to j : 

j 

Pf{O = cl[{-i,i;C\OA;0n,<Pn) 



n=l 



= C'l[un-u{0. (5.2) 

n=l 

Changing the representative from a pair changes the sign of the unspecified numerical 
constants C,C'. The first line of this equation is similar to Majorana's theorem fj3.2|) . 
although in the case of a real function, the decomposition is possible in terms of spin 
1 coherent states, as well as spin 1/2. The second line resembles Maxwell's multipole 
decomposition ()1.2j) . although it is the Majorana function pj(C) that is given here, and it 
is written in terms of the nilpotent vector v{C), not the real vector u{6, 0); furthermore, 
there is no additional summand (this vanishes, since u ■ u = 0). This tranformation will 
be the object of the next section. 
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6. Maxwell's multipoles 



In this section, Maxwell's multipole construction p.2j) for a real function f{6, (p) of spin 
j will be proved based on the decompositions ()5.1|) . ()5.2|) . 

f{0,(p) may be obtained from integrating p/(C) over the ^-plane using the inverse 
kernel K^{e,(l)\C), i.e. 

dCdC* 



;i + icp)^ 



^Y^^ riK ■ KC))KCr ■ ^(^^, 0)), 



(6.1) 



n=l 



where the inverse kernel has been decomposed as the conjugate of ()5.1|) . 

This integral, over all stereographic directions depends on the nilpotent vectors 
f{C)- The C itself is now redundant, and the integral can be taken over all u vectors with 
nilpotence and normalisation conditions provided by 5-functions, written as an integral 
in a form similar to usual coherent state integrals. 



f{OA) = c 



I ud^u* 5{\u\^ - 1)5(1/ • u)5{u* ■ V*) JJ(w„ ■ v){v* ■ u{e, 0)) 

n=l 



(6.2) 



The pair of 5-functions with conjugate arguments are equivalent to pair of 5-functions 
with real independent arguments. 

Now, the overall magnitude of u does not affect the 6, (f) dependence on this integral, 
and the normalisation condition may be relaxed. It is most convenient to replace this 
with a gaussian distribution for the absolute value \u\, giving 



C j dVdV exp(-|i/|V2)(5(i/ • u)6{u* ■ u*) JJ(tx„ ■ 

^ n=l 



u{9,(P)). (6.3) 



This gaussian integral will now be evaluated using standard techniques. 

The pair of 5-functions forcing nilpotence may be replaced by their Fourier 
representation, in the (complex) Fourier variable k. Writing the total integral compactly, 
with V = iv^v*). 



(6.4) 



/(^,0) =C dndn* dVdV*exp(-F*-M- V) JJ(ii„-i/)(i/*-'u(0,0)), 

^ ^ n=l 

where the 6x6 matrix M in the exponent is given by 

/I i/€ \ 

1 iK 

1 1 iK 

4 i/t* 1 

iK* 1 

\ i/t* 1 / 

The inverse M-^ is the same as M with i replaced by — i and divided through by 
(1 + l^p). It is now possible to find the gaussian integral over i/, u*. 



M 



(6.5) 
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This integral is possible by gaussian integration by parts ('Wick's theorem' in 
quantum field theory). Up to a multiplicative constant, the gaussian integral is 
(detM)~^/^ = (1 + times a combinatorial term, which is now described. 

This combinatorial term is a sum over all possible pairings involving u and ly* in the 
product in ()(i.4j) . Each summand is the product of the coefficients from the pairing (in 
terms of the inner product of the vectors tt„ and u{6, cj))) times the component of 
relating the relevant pairing: l/(l + |Kp) for u with i/*, — iK*/(l + for u with and 
-iK/(l + for v* with i/*. The sum is therefore over the product of the dot product 
of the paired vectors, times the appropriate element of M^^ : i.e. w„ ■ u{9^ 0)/(l + l^p), 
iKu{9, (j)) ■ u{9, 0)/(l + I^P), or —iK*Un ■ Un- 

Clearly, there is only one distinct summand which pairs all the u{9, cj)) vectors with 
the Un vectors; there is no k or k* in the numerator of this term. All of the other 
summands have at least one pairing of u{9, cj)) with itself, and two Majorana vectors 
paired. With this in mind, ()6.3|) becomes 

^U{9,<P)-UA j ^^-^^^^^ + U(^,0)-1.(^,0) j dKdK*F{K,K*).{Q.Q) 

F{k, k*) is the sum of all the remaining pairings. The k, k* integrals may now be 
performed; the first summand gives a number, the second, a complicated (but unique) 
function involving u{9,(p) and the u„. This final integral therefore gives ()1.2j) . and the 
Maxwell multipoles have been demonstrated to be exactly the Majorana vectors of the 
function. 

7. Discussion and comparison with alternative proofs 

The main result of this paper is a proof of Sylvester's theorem (|1.2p using Majorana's 
sphere and the algebra of quantum spins. Having started with the real function f{9, (p), 
the strategy was to transform from the fully folded representation (basis functions of 
S0(3)) to the Majorana function (coherent basis functions of SL(2,C)) in the fully 
extended representation, using the extending kernel. Using reality, the Majorana 
factorization is in terms of spin 1 states, rather than the more usual spin 1/2 states; 
for clarity, these were written in vector notation with the spin +1 coherent states 
represented by nilpotent vectors. The transformation back to the folded representation 
was effected by transforming by the inverse kernel, which, as a real state, was also 
factorised. The direction dependence of the ffnal transformation was finally determined 
by transforming the integral on the sphere to a gaussian integral. An additional feature 
of the derivation was the representation of the spherical harmonics via symmetric 
polynomials, between the first and second lines of (jS.lj) . 

The proof provided by Backus |j4j is, in fact, rather similar to the one provided here. 
He constructed the ring homomorphism from the space of homogeneous polynomials of 
order j in three real variables x, y, z, to the space of homogeneous polynomials of order 2j 
in two complex variables ^, rj; the homomorphism he used (|1| equation (50)) is precisely 
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that obtained by using, as cartesian components, quadratic forms in ^ and rj with respect 
to the the Pauh spin matrices times a real antisymmetric matrix. ^ and rj are therefore 
the components of a spin 1/2 spinor, and the mapping is equivalent from going between 
the fully folded and fully extended representations. His factorization of the complex 
polynomial, on replacing ^ — > C ^iid r/ — »• 1 is analogous to Majorana's theorem; reality 
of the original polynomial gives antipodality of the roots. The transformation back 
is made using the first isomorphism theorem for (homogeneous polynomial) rings, by 
mapping into the quotient ring, factoring the original ring in x, y, z with the kernel of 
the spinor mapping, i.e. homogeneous polynomials of the form r ■ rF. In the present 
proof, the transformation is realised concretely by integrating with respect to the inverse 
of the extending kernel. Thus Backus's proof, using spinors in disguise, is equivalent to 
the present one. 

The proof of Zou and Zheng ^ is also similar, but in the language of tensors as 
used in continuum mechanics. Specifically, Sylvester's theorem is applied to traceless 
symmetric tensors (i.e. the unique irreducible tensor operators of maximal order j in the 
decomposition of arbitrary tensors of rank j). As with Backus's and the present proof, 
a complex polynomial is constructed by contracting with a complex tensor function of 
the complex variable (, whose complex roots occur in antipodal pairs. Contracting with 
this tensor function is equivalent to transforming to the fully extended representation. 
A further proof using the same general technique, but the different physical language 
of the spacetime spinor calculus described by Penrose and Rindler |211 ES] , is described 
in [Appendix B In this case, the objects satisfying Sylvester's theorem are symmetric 
hermitian spinors in 3 spatial dimensions (rather than 4 in spacetime). 

Both Courant and Hilbert jH], and Sylvester [2] use identities from the theory of 
algebraic curves, by directly complexifying x, y, z space, and looking for the coincidences 
of the cones f{x,y,z) = with x^ + y"^ + z"^ = 0. Here, Bezout's theorem is used to 
provide the existence of the multipole directions, rather than the fundamental theorem 
of algebra. 

Majorana's sphere construction remains comparatively obscure in the quantum 
mechanics literature, despite providing a natural framework for describing the rotation 
of spins (and, in general, of functions on the sphere). It has been used here, through 
quantum mechanical language, to demonstrate a fundamental property of real states as 
multipoles, which provides a further role for the antipodal pairs of Majorana vectors. 
It is striking to note that this philosophy and property of spinors goes back at least as 
far as Maxwell. 



Acknowledgements 

I am grateful to Gerald Kaiser, Iwo Bialynicki-Birula and Jonathan Robbins for 
discussions, to James Vickers for reading through [Appendix B[ and particularly to John 
Hannay for much insight, and for introducing me to the Majorana sphere. This work 
was supported by the Leverhulme Trust. 



Sylvester's theorem, Maxwell's multipoles and Majorana's sphere 



13 



Appendix A. The equivalence of the two Maxwell multipole representations 

(UlI), (HHD 

It suffices to prove equivalence for a single term, that is, 

^m- = c^+-^, (A.i) 

for p + g + s = j, C some constant independent of p, q and s, and F some other 
homogeneous polynomial in x, y, z of order j — 2. 

To begin, r is replaced by its Fourier integral representation, 

l3j 



dldldll/r = dldldl^^ j i^exp(ifc-r) 



1-^ 



27r2 / A;2 



A;^A:^fc:exp(ifc-r) 



■j poo j- 

— j d\ J d^kkPklKexi){-k^X/2 + ik-r) (A.2) 



47r2 

where in the second line, the partial derivatives have been performed, and in the third, 
itself is represented as an integral over A. The integrals in kx, ky and kz can now be 
evaluated separately; following the appropriate linear transformation k^ = t/V^ + ix/A, 
the kx integral transforms to 

j dkx kPexp{-kl\/2 + ikxx) = '''''^^'fyi^^^ J dt (VAt + ix)''exp(-tV2). (A.3) 

The bracketed term to the pth power in the integrand must now be expanded; all terms 
involving an odd power of t vanish by symmetry, and the nonvanishing terms are real 
due to the i-' factor in the front of ()A.2|) . 

The {ixy summand term may be integrated without problem as a usual gaussian 
(it is independent of p). Multiplying this by the analogous {iy^ and (i^;)* terms, this 
summand of ()A.2|) gives, up to a constant factor 



^V^^ / T4Lexp(-rV2/A) = / dp p^-'/' exp{-p/2) 



= CxPy''zyr^^+\ (A.4) 

In the ffist line, A has been replaced by p = r^/A. This proves the ffist part of ()1.2|) . 

In all of the summands not included in the expression ()A.4|1 . the exponent of A in 
the denominator is less than j + 3/2 by some integer (since only even powers of t in ()A.3|1 
contribute). On transforming to p, the power of r in the resulting denominator is less 
by at least two; l/r^-'"^ therefore factors out, giving the form for the second summand 
in ()1.2j) . Therefore, the two forms of Maxwell's multipole expansion are the same. 



Appendix B. Sylvester's theorem using spacetime spinor calculus 



A further proof of Sylvester's theorem, similar in form to that proved in the main text, 
and lUEj, will be outlined here using the calculus of relativistic spinors, as described by 
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Penrose and Rindler [211 12H] , referred to as I and II in this appendix. It should be noted 
that if algebra of spinors in 3-space were to be used from the outset, as in JH] chapter 
VIII, leads to a more straightforward demonstration (see, in particular, the final section 
of this chapter). 

Using index notation and employing the summation convention, the function f{6, (p) 
with total angular momentum number j may be written 

/(^,0) = r,...x(^,0)---«n^,0), (B.i) 

where u°'{6,(j)) denote the components of the unit vector u{6,(f)) (and, in the following, 
the vector itself). In ()B.1|) . Ta...g = T(^a-g) denotes a real, completely symmetric rank j 
tensor (with round brackets (■ ■ ■) denoting complete symmetrization, as usual), which 
is traceless with respect to contraction with any pair of indices. Where convenient, it 
will be assumed that the basis is fixed. 

It is important to note that to use the spacetime spinor calculus, the indices in (jB.lll 
have to be over 4-dimensional minkowskian spacetime (with signature —2). However, 
Maxwell's multipoles are a feature of 3-dimensional space only, so it is assumed that there 
is an absolute newtonian time direction, that is, a timelike 4-vector v"" = (y^, 0,0,0), 
and Ta...g,u°'{6,(l)) lie on the 3-dimensional spacelike hypersurface orthogonal to f". 

All tensor indices a,b, . . . can be replaced by pairs of spinor indices AA', BB', . . . , 
using the Infeld-van der Waerden symbols g^^'', in the present case, these are simply the 
three Pauli matrices and the identity, with a 2~^/^ prefactor (I section 3.1). The primed 
indices are complex conjugated coordinates, and spinor indices are raised and lowered 
using the antisymmetric symbols eAB,£^^,£A'B',£^'^\ where £01 = ^o'v = = = 
+ 1 (I section 2.5). 

Thus, for instance, 

VaV" = VAA'V^^' = 2, (B.2) 

implying that v^^' is simply the identity matrix \/2g(f'^' . Therefore (II p461), v^\v^, 
may be used to convert between primed and unprimed indices; the relevant components 
are 

v^' = v^, = vl' = vl, = 0, vf = 4 = +1, vl' = v^, = -1. (B.3) 

Since the tensor Ta...g is symmetric and traceless, the corresponding spinor 
Ta---ga'-g' = T(^A-G)(A'-G') is totally symmetric in its primed and its unprimed indices 
(I pl46). Furthermore, since the tensor is real, the spinor is hermitian: 

[T(A-G){A'-G')] = T^A'-G'){A-G) = '^(A-G){A' -G') = T(A-G){A' -G')- (B.4) 

In the second line, the indices have been rearranged (since the ordering of primed indices 
with respect to unprimed is unimportant), and in the third, reality of the components 
of T was used (I ppl23-4). Since T is orthogonal to the timelike vector v, in terms of 
spinors, 

T(A...G){A'...G')^^^' = 0. (B.5) 
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This spinor may be expressed completely in unprimed indices, 

'r(A-G)iAi-Gi) = T(^A-G)(A'-G')vi-, ' ' ' Vq^. (B.6) 

Contracting any two indices from each of the two sets gives 

r{A-G)iAi-Gi)£ - ^(A-G)(A'-G')^Ai^ ^Bi " ' " ^Gi 
- -J- {A-G){A'-G>)V %i ■ ■ ■ %i 

= (B.7) 
by ()B.5j) . This shows that r is symmetric with respect to indices in its two symmetrised 
sets, and thus is totally symmetric: T(^a-g){a^-g^) = ^(a-gAi-Gi)- 

This totally symmetric spinor may be canonically decomposed (I proposition 
(3.5.18)), 

TiA-GA^-Gi) =a(A---lGa'A,---lG,)^ (B-8) 

where spinors are primed when their index has subscript 1. This step is, of course, 
equivalent to Majorana's theorem, which is proved (I pl62) by using the fundamental 
theorem of algebra on the order 2j (Majorana) polynomial 

r{C) = nA-G,)^ ■■■e^ where = J ) • (^-9) 

and a spin frame is chosen where the complex number C 7^ 0. 

It is necessary to show that the roots of this polynomial are antipodal, since then 
the spinors in the decomposition ()B.8|) are arranged in antipodal pairs, which will be 
demonstrated similarly to ()3.6p . The spinors may be replaced by conjugated, primed 
spinors t)*"^' using 

V*^' = -v^^^- (B.IO) 
From defined in (p^ . (pin|l gives 

The polynomial ()B.9jl may thus be written in terms of ^ and rj, 

r(C) = (-l)^V..G)(A'...G')e^---eV^'---r/*^'. (B.12) 
Since the components of T are real, conjugation of the polynomial gives the same result 
as in (jH.fjj) : 

r{cr = {-iyTiA'...G')iA...G)C^' ■ ■ ■ r^'^^ ■ ■ ■ 

= {-iynA-G)iA'-a)V''---ri''C'''---C''' 

ocr(-l/C). (B.13) 

Thus the roots are antipodal. 

Transforming the Maxwell multipole form of ()1.2|) into spinor form is equivalent 
to Ti^A...G^)i in the sense of I pl40: it equals it, plus combinations of lower rank spinors 
combined with the antisymmetric symbol. This equivalence relation is equivalent to 
ignoring the second summand in the multipole decomposition p.2|) . Thus Sylvester's 
theorem is proved. 
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